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ZJ ' Abstract. We review the theory and observations of star cluster disruption. 

, The three main phases and corresponding typical timescales of cluster disruption 

are: I) Infant Mortality (~ 10 7 yr), II) Stellar Evolution (~ 10 8 yr) and III) 
C*") | Tidal relaxation (~ 10 9 yr). During all three phases there are additional tidal 

external perturbations from the host galaxy. In this review we focus on the 
physics and observations of Phase I and on population studies of Phases II 
& III and external perturbations (concentrating on cluster-GMC interactions). 
Particular attention is given to the successes and short-comings of the Lamers 
cluster disruption law, which has recently been shown to stand on a firm physical 
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1. Introduction 

The vast majority (perhaps all) of stars are formed in a clustered fashion. How- 
ever, only a very small percentage of older stars are found in bound clusters. 
These two observations highlight the importance of clusters in the star-formation 



process and the significance of cluster disruption. The process of cluster disrup- 
tion b egins soon after, or concurrent with, cluster formation. lLada fc Ladal 
(|2003h found that < 10 % of stars formed in e mbedded c l usters end up in bound 
clusters after ~ 10 8 yr. IWhitmord (|200.t ) and lFall et aD OoO-V) have shown that 



at least 20%, but perhaps all, star formation in the merging Antennae galaxies 
is taking place in clusters, the majority of which are likely to become unbound. 
The case is similar in M51, with > 60% of all young (< 1 Myr) clusters likely 



to be destroyed withi n the first 10s o f Myr of their lives ([Bastian et al.l 12005). 
On longer timescales, lOortJ (ll95ST ) and lWielenl (|l97lT ) noted a clear lack of older 
(> few Gyr) open clusters in the solar neighbourhood and lBoutloukos &: Lamersl 
( 2003) found a strong absence of older clusters in M51, M33, SMC, and the solar 
neighbourhood. 

The lack of old open clusters in the solar neighbourhood is even more 
striking when compared with the LMC, which contains a significa nt number 
of 'bl ue globular cluster s' with ages well in excess of a Gyr (e.g. iGascoignel 
119661 : Ide Griis fc Andersl l2006h . This difference can be understood either 
difference in the formation history of clusters or as a differ e nce in the disruption 
timescales. This later scenario was suggested bv iHodgel (|l987l ). who directly 
compared the age distribution of Galactic open clusters and the SMC cluster 
population. He noted that there are 10 — 15 times more clusters with an age of 
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1 Gyr in the SMC as compared to the solar neighbourhood (when normalising 
both populations to an age of 10 8 yr) and concluded that disruption mechanisms 
must be less efficient in the SMC. 

Much theoretical work has gone into the later scenario, with both ana- 
lytic and numerical models of cluster evolution predicting a strong influence of 
the galactic tidal fiel d on the dissolution of star clusters (for a recent review 
see IBaumgardtl 12006"), Only recently has there been a large push to under- 
stand cluster disruption from an observational st andpoint in various external 
potentials, making e xp licit comparis on wi th models dBoutloukos Lamersl2003l : 
lLamers et aLll2005al lbl: iGieles et al.ll2005l : lLamers fc Gielesl 1200^1 . " 

We direct the reader to the review by Larsen in these proceedings for a 
historical look at the observations and theory of cluster disruption. 

1.1. Phases of cluster disruption 

While cluster disruption is a gradual process with several different disruptive 
agents at work simultaneously, one can distinguish three general phases of clus- 
ter mass loss and disruption. As we will see, a large fraction of clusters gets 
destroyed during the primary phase. The main phases and corresponding typ- 
ical timescales of cluster disruption are: I) Infant Mortality (~ 10 7 yr), II) 
Stellar Evolution (~ 10 8 yr) and III) Tidal relaxation (~ 10 9 yr). During all 
three phases there are additional tidal external perturbations from e.g. giant 
molecular clouds (GMCs), the galactic disc and spiral arms that heat the cluster 
and speed up the process of disruption. However, these external perturbations 
operate on longer timescales for cluster populations and so are most important 
in Phase III. In Fig. ^ we schematically illustrate the three Phases of disrup- 
tion and the involved time-scales. Note that the number of disruptive agents 
decreases in time. 

In this review we will focus on the physics and observations of Phase I as 
well as on recent population studies aimed at understanding Phases II and III 
on a statistical basis. For a recen t revie w on the physics of Phases II and III, 
we refer the reader to IBaumgardtl ( 20061 ) . 



Before proc eeding, it is worthwhile to consider our definition of a cluster. 
ISchweizerl 0006) defines a cluster to be a gravitationally bound stellar association 



which will survive for 10-20 crossing times. This definition implies that the stars 
provide enough gravitational potential to bind the cluster and ignores the role 
of gas in the early evolution of clusters. In this review, we will define a cluster 
as a collection of gas and stars which was initially gravitationally bound. The 
reason for this definition will become evident in Section |2~1 



2. Infant Mortality 

Rece nt studies on the populations of young star clusters in M51 dBastian et all 
2005) and the Antennae galaxies (|Whitmorell2003l : iFall et al.l [2005l have shown 
a large excess of star clusters with ages less than ~10 Myr with respect to what 
would be expected assuming a constant cluster formation rate. The fact that 
open clusters in the solar neighbourhood display a similar trend ( Lada &: Ladal 
2003) has led to the conclusion that this is a physical effect and not simply 
that we are observing these galaxies at a special time in their star-formation 
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Disruption phases 

I) Infant mortality Cas expulsion 
► 

II) Stellar evolution Mass loss due to stellar winds 
► 

III) Tidal relaxation Two body relaxation 
► 

External perturbation Tidal heating (GMCs, disc crossing, etc.) 



io u 10 io u icr 10 

Time scale (yr) 

Figure 1. Schematic overview of the three phases of cluster disruption con- 
sidered and the responsible physics that drives the disruption. 



history. If one adopts this view, then we are forced to conclude that the majority 
(between 60-90%) of star clusters become unbound when the remaining gas (i.e. 
gas that is left-over from the star formation process) is expelled. These clusters 
will survive less than a few crossing times. 

2.1. Gas expulsion 

Suppose that a star cluster is formed out of a sphere of gas with an efficiency e, 
where e = M s t ar s/(^ stars + M gas ). Further suppose that the gas and stars are 
initially in virial equilibrium. If we define the virial parameter as Q = —2T/W, 
with T the kinetic energy and W the potential energy, virial equilibrium implies 
Q = 1. Finally, suppose that the remaining gas is removed on a timescale faster 
than the crossing time of stars in the cluster. 

In such a scenario the cluster is left in a super-virial state after the gas 
removal, with Q = 1/e, and the star cluster will expand since the binding energy 
is too low for the stellar velocities. The expanding cluster will reach virial 
equilibrium after a few crossing times, but only after a (possibly large) fraction 
of the stars have escaped. This process has been shown to remove a significant 
amount of the stellar mass of a cluster, and if e < 0.3 the entire cluster w ill 
become unbound on a timescale of 10s of Myr (lTitiikovlll973:lGoodwinlll997al lb1: 
iKrouna fc Boihill2002t iBoilv fc Kroi7pl]l2003al lbl: iBastian fc Goodwinll20o'rf L 

Rapid gas removal of the type discussed above leaves distinct observables. 
In Figure |21 we show the surface brightness profiles of three young clusters (left 
panels) as well as two results of iV-body simulations (right panels) of clusters 
including the effects of rapid gas removal. All three young clusters show an 
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exces s of light at lar ge radii with respect to the best fitting EFF ( Elson et al.l 



11987ft or iKind (|l962T l profiles. This is in good agreement with the predictions 



of the simulations, in which an unbound halo of stars is removed (although still 
appearing to be associated with the cluster for 10s of Myr) du e to the rapid 
change of the gravitational potential ( Bastian k, Goodwinl 120 06). Such excess 



light at large radii has al so been found in young clust ers in the Antennae galaxies 
(Whit more et al.lll999l ). iGoodwin fc Bastianl (2006) show that for values of e of 



0.1 and 0.6, clusters will lose 75% and 10% of the stellar mass respectively within 
the first ~ 20 Myr of their lives. 

Thus we see that this is an extremely efficient way to rapidly disperse 
stars from young clusters into the field. This mechanism provides a natural 
explanation for the observed diffuse UV li ght in the field of starburst galaxies 
( Tremonti et~aTll200ll : iGhandar et al.ll2005T ) and supports the scenario of these 



authors that this light is due to rapidly dispersing young clusters. 

Whether or not a cluster survives this phase, and hence more than 10-20 
crossing times, is largely dependent on the star-formation efficiency of the GMC 
core in which the cluster formed. Thus, two clusters with exactly the same pa- 
rameters (radius, mass, metallicity, external potential field, etc) may experience 
two radic ally different evolutionary p aths if their star-formation efficiencies are 
different. IGoodwin Sz Bastianl ((20061 ) have used the internal dynamical proper- 
ties of young clusters in order to estimate their e- values. No clear trend of e on 
cluster (stellar) mass or radius was found. 

2.2. Interpretation of the properties of young clusters 

Even if a cluster survives the gas removal phase, this phase can significantly effect 
the observed properties of the cluster. Henc e, deducing the initial p roperties of 



tne observed properties ol tne cluster. Hence, deducing tne initial properties ol 
a cluster from its current state is not trivial. iKrouoa Sz Boilvl ( 20021 ) have noted 



the strong effe ct of residual gas removal o n inferring the initial stellar mass of 
a cluster, while IGoodwin Bastianl (|2006l ) have refined the mass loss estimates 
and shown that measurements of the current radii of young clusters may not 
reflect the initial nor the final value. Additionally. IGoodwin Bastianl (2006) 
show that this effect can mimic stellar IMF variations in young clusters. 



3. Population Studies of Cluster Disruption 

The clusters that have survived the gas removal phase are subject to disruption 
Phases II and III (§ ll.l.j) as well as tidal effects. Disruption due to these effects 
can be studied on individual clusters, of which the recent obser vations of the 
dissolving globular clusters Palomar 5 ( Qdenkirchen et al.ll200il ) are probably 



the most spectacular example. However, much can be learned by approaching 
this problem from a cluster population point of view. 

3.1. The Lamers disruption law 

Suppose that clusters are formed continuously with a constant cluster formation 
rate (a constraint which we can relax later). Also, we will assume that we 
know the cluster initial mass function (here taken to be a p ower law of the form 
N d dM c i oc M c r a dM c i with a = 2 fe.g. lZhang fc FaJJIl999t Ide Griis et all 2003) 
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Figure 2. Surface brightness profiles for three young clusters (left - M82-F, 
NGC 1569-A, and NGC 1705-1) and two TV-body simulations which include 
the rapid removal of gas which was left over from a non-100% star- formation 
efficie ncy (right). The s olid (red) and dashed (blue) lines are the best fitting 
EFF ijElson et alJll987j) and King l)Kindll962|) profiles respectively. Note the 
excess of light at large radii with respect to the best fitting EFF profile in 
both the observations and models. This excess light is due to an unbound 
expanding halo of stars caused by the rapid ejection of the remaining gas after 
the cluster forms. Hence, excess light at large radii strongly implies that these 
clusters are not in dynamical equilibrium. For details of the modeling an d 
observations see lBastian fc Goodwin! l|200fiD and lGoodwin fc Bastianl ((2006) . 



and that clusters can be detected down to a known magnitude limit. Finally, we 
will assume that the disruption time of a cluster depends on the cluster mass, 
such that more massive clusters survive longer (on average) than lower mass 
clusters. For this final assumption we will adopt a function of the form: 

i dis = t 4 (M cl /10 4 M o )^, (1) 
wher e t± is the disruption time of a 10 4 M cluster and < 7 < 1 ( Boutloukos Lamersl 



2003). The beauty of this formulation is that it only has two variables, £4 and 
7, and as we will see, provides extremely good fits to observations. 

3.2. Application to various cluster populations 

The formulation provided above, when combined with the given assumptions, 
allows for the parameters £4 and 7 to be found from age and mass distribu- 
tions of clusters. The first survey using this technique was carried out by 
iBoutloukos k, Lamersl ( 20031 ) on cluster populations in M51, M33, the SMC, and 
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Figure 3. Comparison between the mass loss following from the TV-body 
simulations of clusters with different number of stars, different concentration 
and on different orbits (left). The mass loss due t o stellar evolu t ion is n ot 
shown. In the right panel the analytical model of lLamers et al.l l|2005al) is 
shown. 



the solar neighbourhood. They made a sudden disruption assumption, meaning 
that the cluster is in the sample with its initial mass until idis; when it is dis- 
rupted. The somewhat surprising result from this study was that, while 7 had 
more or less the same value in all environments studied ((7} = 0.62 db 0.06), £4 
varied by over two orders of magnitude, with values of ~ 80 Myr in the central 
regions of M51 to ~ 8 Gyr in the SMC. 

The simple sudden disrupt ion assumption was improved in a more recent 
model by lLamers et al. I (|2005ari . where a gradual loss of cluster mass was im- 
plemented. They assumed that the cluster mass decreases exponentially with a 
time-scale that decreases as the cluster mass decreases. This is done by saying 
that the mass loss per unit time (M) relates to idis as: 

M = Md/tdta oc M cl x ~\ (2) 

with idis from Eq. ^ This very simple analytical description for cluster mass 
loss shows remarkably good agreement when compared to the mass los s following 
from the detailed iV-body simulations of lBaumgardt fc Makind (|2003T t . In Fig. 



we show a direct comparison of M from the iV-body simulations of clusters with 
different density pr ofiles and on different orbits (left) and the above mentioned 
analytical model of lLamers et al.l ( 2005al ) (right). 



In both graphs the time is normalised to tdis and only the mass loss due to 
stars escaping the cluster is shown, i.e. mass loss due to stellar evolution (SEV) 
is not shown. In addition, there is a coupling between the two types of mass 
loss: if stars loose mass, the cluster wil l expand and more stars are pushed over 
the tidal boundary. The simulations of lBaumgardt fc Makinol l|200.^ considered 



SEV, therefore, their M results shown in Fig. |3] do include tidal M induced by 
SEV. For this reason we can simply add the mass loss due to SEV, taken from 
an SSP model, to Eq. EJ 

In a series of follow-up works, it has been shown that the similarity of the 
value of 7 in various environments strongly implies a uniformity in the cluster 
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disruption process, while the varying values of t$ is due to the different tidal field 
strengths (and gas contents) of the galaxies studied. Galaxies with strong tidal 
fields, as, fo r example, derived fr om their rotation curves, having shorter disrup- 
tion times (|Lamers et al.l l2005bt) . Comparison with results of realistic iV-body 
models performed bv iBaumgardt Makind (l2003h have placed this empirical 
disrupt ion law on a solid physical footing (|Lamers et alJl2005al ). lLamers et al.l 
(2005a) have also derived a formula for the predicted mass and age distributions 
of cluster samples that includes both stellar evolution and disruption for any 
star formation history. 

iGieles et al.l (|2005h inserted the Lamers disruption law into a cluster popu- 
lation synthesis model. This method has two distinct advantages over the earlier 
formulations. The first is that it removes the requirement of a constant cluster 
formation rate, and second, it uses the age and mass distributions together to 
find 7 and £4. The case of M51 is shown in Fig. One first begins by construct- 
ing an observed number density grid in age-mass space (upper-left panel where 
the shading corresponds to the logarithm of the number of clusters found within 
that cell). Then one generates a large number of models with different values of 
£4, 7, (time dependent) cluster formation rates, etc. and compares these models 
with the observed grid. The resultant xt diagram is shown in the bottom panel 
of Fig. The best fit model is shown in the top right panel of Fig. |1J 

T his cluster population synthe sis (CPS) technique, also used in a similar 
way by Dolp hin Sz Kennicuttl (|2002f ) to derive the properties of the cluster popu- 
lation in the galaxy NGC 3627, holds great promise in disentangling the myriad 
of effects present in cluster populations. In principle, the dependences of cluster 
size, galactocentric radius, star-formation efficiency dependent infant mortality 
rates, or alternative cluster disruption models can be taken into account by this 
technique. For this technique to be fully exploited one needs large samples of 
cluster populations with known ages and masses. Datasets suitable for these 
kind of studies are beginning to be collected and released. Several face-on spiral 
galaxies have been ima ged in multiple filter s with the hig h resolution /wide fi eld 
HST/A CS camera fe.g lBarmbv et alJl200fil for M101 and lGieles et~atll2006bl for 
M51). 



3.3. Possible objections 

It is worth noting possible objections to the Lamers disruption law. The first 
comes from lFall et all {2005) who find that in the Antennae galaxies the number 
of clusters decreases in time (r) as dN/dr cx r , independent of cluster mass. 
This may be explained if the disruption timescale £4 due to tidal field effects 
(e.g. Phase II & III) is greater than or similar to the maximum age in the 
sampl e. The cluster d isruption due to tidal effects would not yet be present 
m the Fall et alJ (|2005l ) sample, instead the decrease in cluster numbers would 
be the result of infant mortality and the fading of clusters. St udies of infant 
mort ality in M51 also suggest that the effect is mass independent (fBastian et al.l 
2005). In fact, if infant mortality was not (mostly) independent of cluster mass 
we would expect the embedded cluster mass function to be significantly different 
from the optically selected cluster mass function. 

In Fig. [21 we show the dependence of the mass function slope of a multiple 
age cluster population on the ratio of the t± and the maximum age of the cluster 
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Figure 4. Top left panel: The observed age-mass diagram of clusters in 
M51. The shading represents the logarithm of the number of clusters found 
within the corresponding box in age-mass space. The solid line represents the 
detection limit of the sample. Top right panel: The best fitting cluster pop- 
ulation synthesis model which includes the detection limit of the observations 
and has fit on 7 and t±. Bottom: The xl diagram in 7-^4 space. The best 
fitting model is marked with an 'x', while the accepted fits (xl < xl best + 1) 
are shaded. Reproduced from lGieles et al.l f|200f5f) . 



in the sample (t m ax)- Clusters were created continuously over 1 Gyr with an 
initial mass function of a power-la w with index —2. T he Lamers disruption law 
was applied in the same way as in lGieles et alJ ((2005). The important thing to 
take away from this figure is that as t<jis approaches and exceeds i max the mass 
function is less affected by disruption and s o it retains its in itial form, i.e. the 
right panel in Fig. Improbably applies to the lFall et alJ (|2005T ) sample. 

A second observa t ion se emingly contradicting the Lamers disruption law is 
that o f lChandar et"aL1 {2006) who find an interm ediate age (~ 4—7 Gyr) globular 
cluster in M33. In M33. lLamers et alJ (|2005bft find a t A value of ~ 600 Myr, 
implying a disruption time of ~ 2.5 Gyr for a 10 5 Mq cluster . However, as the 
authors note, the value of £4 derived bv lLamers et alJ (|2005bf l was presumably 
of the thin disk of the galaxy, and if the intermediate-age cluster is part of 
the thick-disk or halo of the galaxy then the expected value of £4 would be 
significantly larger t han that quoted . Add itionally, it should also be noted that 
the mass derived by IChandar et alJ (|2006h is the present mass of the cluster. 
The cluster presumably started with a much higher mass and disruptive effects 
have brought this cluster into its current state. If the present mass of the cluster 
is 1 x 10 5 Mq, then its initial mass (after infant weight loss) would have been 
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Figure 5. Representation of how the slope of the mass function changes 
as a function of the ratio between the disruption time (t^ is ) and the age of 
the oldest cluster in the sample (in this case assuming a constant cluster 
formation rate) . Note that if idis approaches or exceeds t max the slope of the 
mass function approaches the initial mass function value (i.e. —2). 
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Figure 6. Comparison between th e observed (filled circles) and predicted 
(solid line; from the TV-body models of lBaumgardt fc Makinol20f)3t) disruption 
time of a 10 4 Mq cluster, £4, as a func tion of the mean densit y p a mb hi M©pc~ 3 
of the host galaxy. Reproduced from lLamers et all l)2005b|) . 



5 x 10 5 M.(7) (assuming an age of 5 Gyr), using the value of £4 for M33 given by 



b x ilr Mm ( a ssuming 
lLamers et all (|2005bf) . 



4. External Disruption Effects: GMCs and Spiral Arms 

4.1. The disruption time due to external perturbations 

As more and more galaxies (and environments) have their characteristic disrup- 
tion timescales measured, it is useful to compare the results to iV-b ody models in 
order t o check for consistency between the two. This was done in lLamers et al.l 
(2005b) who compared the t± values derived for the SMC, M33, M51 and the 
solar neighbourhood to the iV-body models of lBaumgardt fc Making 1)20031 ) and 
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Figure 7. Snapshots of a star cluster undergoing a head-on encounter with 
a GMC of M„ = 10 4 M Q with V max = 20cr rms . Top: The motion of the 
GMC is along the x-axis and the line of sight is perpendicular. Bottom: The 
motion of the GMC is into the page and the line of sight is along the GMC 
trajectory. The arrows in the left-hand lower corner of the left-hand panels 
are parallel to the direction of motion of the GMC. The GMC is shown with 
grey shades based on the surface density of a GMC with a„ = 5.8a c . The 
time with respect to the moment of e ncounter is indicated in each panel of 
the top row. See iGieles et alJ l)2006cj) for a description of the methods and 
parameters used. 



Portegies Z wart et al. ( 19981. 12002ft which sample a large range in the ambient 
densities of the host galaxies. Their results, shown in Fig. El are intriguing. 
While the predicted and observed disruption time of the SMC are in excellent 
agreement, the disruption times of the Galaxy, M33 and M51 are observed to 
be much shorter than predicted by iV-body models. This result is particularly 
surprising given the fact that the mass loss predictions of a single cluster are 
in excellent agreeme nt between the Lame rs empirical description and that given 
by iV-body models ( Lamers et al.ll2005al and Fig. |3J). 

Thus, we are left to ask, what physical effects are not included in the N- 
body models that may be responsible for disrupting clusters? The iV-body 
models used in the comparison were carried out in a smooth logarithmic potential 
which does not reali stically represent the thin disk components of disk galaxies. 
IGieles et al.l (|2Q06al O) have attempted to add encounters with giant molecular 
clouds (GMCs) and spiral arm passages to the iV-b qdy models. In Fig . [7] we 
show an example of a cluster-GMC encounter ( from IGieles et all 120060 ). The 
parameters of this run are for typical open clusters and GMCs in the solar 
neighbourhood. The top panels show an edge-on view for five different time 
steps, while the bottom panels show a view along the trajectory of the GMC. 

Encounters with GMCs present the most important ext ernal perturbation 
which cause mass loss of star in clusters. IGieles et al.l (20063) find that idis due 
to encounters with GMCs scales as 
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M cl \ /3.75pc x 3 



where 5 = 1 for the solar neighbourhood and 5 scales with the surface density 
of individual GMCs (S n ) and the global GMC density (p n ) as S oc (S n p n ) _1 . 
The scaling of S with p n implies that it does not matter if the molecular gas 
is distributed over a large number of low mass clouds or a small number of 
massive (giant) clouds. This makes it easy to estimate t^is from the observed 
molecular gas density. Indeed, for M51, where the molecular gas density is about 
an order of magnitude higher than in the solar neighbourhood, a i^is from Eq. |3] 
of 150 Myr is predicted. This cor responds well with the value derived from 
observations of t 4 = 100 - 200 Myr (|Gieles et aUbOOffl . 



Note that Eq. |3] implies a scaling of t^ s with the cluster density (p c i). This 
seems different than the dependence with M c \ discussed before. However, there 
is only a very shallow relation o bserved betw een cluster half-mass radius (r^) 
and M ch of the form r h oc M cl 01 (|Larsenll20"olh . With this relation, and Eq.[3 it 
follows that for external perturbations t^is Mcl ' 7 \ i- e - very close to the index 
of 7 ~ 0.6 found from observations discussed in § 13.2.1 This suggests that the 
disruptive effect of the tidal field and additional external perturbations can be 
added linearly, resulting in a tdis that depends on M c \ as tdis cx M c i 0,6 . This can 
explain the large variation fou nd in the t&, value derived fr om observations and 
the almost constant 7 = 0.6 ( Boutloukos Sz Lamersl 120031 ) . In §E]we discuss 



some of the pitfalls of these results. 

4.2. Application to the open clusters in the Galaxy 

As seen in the proceeding sections, the observed disruption time of star clusters 
in the solar neighbourhood is a factor of ~ 5 shorter than predicted by iV-body 
models. The inclusion of spiral arm passages and GMC encounters into N- 
body models is a promising way to bring t he predictions into agre ement with 
the observations. This was recently done bv lLame~ fc Gielesl (|2006fl who found 
excellent agreement after the inclusion GMC encounters and spiral arm passages. 
They assume that the different mass loss effects (stellar evolution, tidal field and 
external perturbations) can be added lin e arly. Using t he mass-radius relat ion of 
§ 14.1.1 and the results from iGieles et alJ (|2006al ) and iGieles et alJ (|2006d l they 
analytically model the mass loss due to di fferent effects analytica lly. This is 
illustrated in the left panel of Fig. |S] (from lLamers fc Gielesl 12006^ / Based on 
this mass loss description, the age distribution of open clusters in the solar 
neighbourhood can be predicted (instead of fitted, as was done hitherto). The 
results are shown in the right panel of Fig. |HJ 



5. Discussion 



We showed in I3~&l4~l that the simple Lamers disruption law can successfully 
explain the age and mass distribution of young star clusters populations. Here 
we will discuss other observations lending support to the Lamers law and some 
of the standing problems and uncertainties of this scenario which need further 
attention. 
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Figure 8. Left: The mass evolution of a 10 4 Mq cluster due to various 
disruptive effects. Right: Comparison between t he observed age distribu- 
tion of open clusters (fr om lKharchenko et al.ll2005l) and the predictions from 
lLamers fc Gielesl l|200f)|) for three different maximum masses. 



5.1. Independent checks on the disruption law 

It is reassuring to see that different datasets of various cluster populations all 
com e to similar conclus i ons re garding the disruption of clusters. 

Ide Griis fc Andersl (|200fl) use a variety of studies to look at the cluster 
population of the LMC. They also find a lack of old clusters (with respect to 
what would be expected from a continuous cluster formati on rate) and derive 7 = 
0.56, again in agreement with ot her galaxies studied by iBoutloukos &: Lamerd 
(2003]) and lLamers et all (|2005al ). Note that a lower value of 7 is expected to 
be observed when the typical tdis is of the same order as the oldest clusters in 
the sample (Fig. [SJ), as is the case in LMC. 

Outside th e local group, t he st rongly interacting galaxy NGC 6745 has 
been studied bv Ide Griis et alJ (|2003l ) who found evidence for mass dependent 
disruption, with 7 = 0.54. The rich cluster system of the intermediate-age 
merger remnant NGC 1316 shows a clear bimodal colour distr ibution, with the 
red co mponent presumably being formed during the merger. iGoudfrooii et all 
(2004) showed, using deep HST-ACS images that if one breaks the red compo- 
nent into 'inner' and 'outer' regions (with respect to the galactic centre), that the 
outer region is a continuous power-law while the inner region shows a power-law 
behavior at the high luminosity end and a flattening at the low luminosity end. 
The authors interpret this as evidence for mass-dependent cluster disruption, 
although no attempt was made to find the characteristic disruption timescale or 
the value of 7. 

One standing problem with the Lamers disruption law, also present in other 
studies on disruption, is whether or not an initial power-law cluster initial mass 
function (CIMF) can be transformed into a log-normal distribution, which is 
observed for old globular cluster populations. The Lamers law can create such 
a turnover, however the precise value of the tu rnover mass should be dependent 
on the ambient density ((Lamers et al.ll200 5b). meaning that cluster disruption 
should be more efficient in the inner regions of a galaxy than in the outer regions. 
Thus, without fine tuning the models (e.g. having the same disruption time at 
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all radii due to large radially dependent velocity anisotropies) one would expect a 
radially dependent turnover peak in the globular cluster mass function, which is 
not observed. For a more detailed description of this pr oblem, see the revie w by 
Larsen in these proceedings. Additionally, as noted bv I Waters et al. (200$ the 
Lamers disruption over-predicts the number of low-mass clusters when applied 
to old globular cluster populations. 

5.2. Caveats in the theoretical underpinning 

In § !4.1.l we showed that the scaling of t<jis with M c \ is a power-law with exponent 
7 ~ 0.6. This scaling is similar for two-body evaporation in a tidal field with 
external perturbations, such as shocks by GMCs and spiral arms, and agrees 
well with the observations. However, there are still some caveats in the theory 
explaining this, mostly coming from questions regarding the initial conditions of 
the simulations. 



The first caveat stems from the relation between initial mass and radius 
of the clusters u sed in the simulations. If we p arameterize this relation as 
rh oc M c i A , then iBaumgardt Makinol ( 2003h use A = 1/3, implying that 



their clusters fill their tidal ra dius. However, observations impl y that A = 
— 0.1 (with a large scatter) (|LarserJl2004l : iBastian et al.ll2005l h implying 
that is mostly independent of mass. This shallow relation implies that 
massive clusters are not filli ng their tidal radius, which w ould change the 
dependence of t<jis with M c \ (jTanikawa Fukushigdl20C)5h . 

2. In the derivation of 7 for external shocks (§ 14.1. j) . only clusters in isolation 
were considered. How would the presence of a tidal field affect this result? 

3. How does the relation for t^is change if there exists a relation between the 
concentration parameter and mass of a cluster (i.e. as seen in the Galactic 
GCs reported by Larsen in these proceedings)? 

4. Could an initial mass-radius relation with A = 1/3 be erased during the 
gas removal phase? 

5. What is the effect of the initial mass/luminosity profile used in the simu- 
lations and how does it evolve? e.g. are clusters born with EFF profiles 
which are converted into King profiles? Does the cluster concentration 
alter its mass loss evolution? 

6. How do the external perturbations and the galactic tidal field cooperate? 
Can the mass loss due to both effects simply be added linearly? 



Acknowledgments. Both NB and MG are grateful to Henny for all of his 
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References 

Barmby, P., Kuntz, K. D., Huchra, J. P., & Brodie, J. P. 2006, AJ, 132, 883 
Bastian, N. & Goodwin, S.P. 2006, MNRAS, 369, L9 

Bastian, N., Gieles, M., Lamers, H.J.G.L.M., Scheepmaker, R. A., & de Grijs, R. 2005, 
A&A 431, 905 



14 



Nate Bastian and Mark Gieles 



Baumgardt, H., & Makino, J. 2003, MNRAS, 340, 227 

Baumgardt, H. 2006: Globular Clusters: Guides to Galaxies (astro-ph/0605125) 
Boily, CM. & Kroupa, P. 2003a, MNRAS, 338, 665 
Boily, CM. & Kroupa, P. 2003b, MNRAS, 338, 673 
Boutloukos, S.G. & Lamers, H.J.G.L.M. 2003, MNRAS, 338, 717 

Chandar, R., Leitherer, C, Tremonti, C.A., Calzetti, D., Aloisi, A., Meurer, G.R., & 

de Mello, D. 2005, ApJ, 628, 210 
Chand ar, R., Puzia, T.H., Sarajedini, A., Goudfrooij, P. 2006, ApJL, in press 

(jastro-ph/06064ll 

de Grijs, R., Anders, P., Bastian, N., Lynds, R., Lamers, H.J.G.L.M., O'Neil, E.J. 2003, 

MNRAS, 343, 1285 
de Grijs, R., & Anders, P. 2006, MNRAS, 366, 295 
Dolphin, A.E., & Kennicutt, R.C 2002, AJ, 124, 158 
Elson, R.A.W., Fall, M.S., & Freeman, K.C 1987, ApJ 323, 54 (EFF) 
Fall, S. M. & Zhang, Q. 2001, ApJ, 561, 751 
Fall, S.M., Chandar, R. & Whitmore, B.C. 2005, ApJ, 631, L133 
Gascoigne, S. C. B. 1966, MNRAS, 134, 59 

Gieles, M., Bastian, N., Lamers, H.J.G.L.M., & Mout, J.N. 2005, A&A, 441, 949 
Gieles, M., Athanassoula E., Portegies Zwart S. F., 2006b, MNRAS, submitted 
Gieles, M., Larsen, S. S., Scheepmaker, R. A., et al. 2006b, A&A, 446, L9 
Gieles, M., Portegies Zwart, S.F., Baumgardt, H. Athanassoula, E., Lamers, 

H.J.G.L.M., Sipior, M., Leenaarts, J. 2006c, MNRAS, 371, 793 
Goudfrooij, P., Gilmore, D., Whitmore, B.C., & Schweizer, F. 2004, ApJ, 613, L121 
Goodwin, S.P. 1997a, MNRAS, 284, 785 

Goodwin, S.P. 1997b, MNRAS, 286, 669 

Goodwn, S.P. & Bastian, N. 2006, MNRAS, in press (astro-ph/0609477) 
Hodge, P. 1987, PASP, 99, 724 

Kharchenko, N. V., Piskunov, A. E., Roser, S., Schilbach, E., & Scholz, R.-D. 2005, 

A&A, 438, 1163 
King, I. 1962, AJ 67, 471 

Kroupa, P. & Boily, CM. 2002, MNRAS, 336, 1188 
Lada, C.J. & Lada, E.A. 2003, ARA&A, 41, 57 

Lamers, H.J.G.L.M., Gieles, M., Bastian, N., Baumgardt, H., & Kharchenko, 

N.V. 2005a, A&A, 117, 129 
Lamers, H.J.G.L.M., Gieles, M., & Protegies Zwart, S.F. 2005b, A&A, 429, 173 
Lamers, H.J.G.L.M. & Gieles, M. 2006, A&A, 455, L17 
Larsen, S. S. 2004, A&A, 416, 537 
Odenkirchen, M., et al. 2001, ApJ, 548, L165 
Oort, J. H. 1958, Ricerche Astronomiche, 5, 507 
Portegies Zwart, S.F., Hut, P., & Makino, J. 1998, A&A, 337, 363 
Portegies Zwart, S.F., Makino, J., McMillan, S.L.W., & Hut, P. 2002, ApJ, 565, 265 
Schweizer, F. 2006: Globular Clusters: Guides to Galaxies (astro-ph/0606036) 
Spitzer, L. 1987, Dynamical evolution of globular clusters (Princeton, NJ, Princeton 

University Press, 1987, 191 p.) 
Tanikawa, A., & Fukushige, T. 2005, PASJ, 57, 155 

Tremonti, C.A., Calzetti, D., Leitherer, C, & Heckman, T.M. 2001, ApJ, 555, 322 
Tutukov, A.V. 1978, A&A, 70, 57 
Vesperini, E. 1998, MNRAS, 299, 1019 
Vesperini, E. & Zepf, S. E. 2003, ApJ, 587, L97 

Waters C.Z., Zepf, S.E., Lauer, T.R, Baltz, E.A., & Silk, J. 2006, ApJ, in press 

lastro-ph/0607231 
Wielen, R. 1971, A&A, 13, 309 

Whitmore, B.C. Zhang, Q., Leitherer, C, et al. 1999b, AJ, 118, 1551 
Whitmore, B. C. 2003, in A Decade of HST Science, eds. Mario Livio, Keith Noll 
& Massimo Stiavelli, (Cambridge: Cambridge University Press), 153. (widely 



Theory and Observations of Cluster disrupts 



referenced as Whitmore 2000, astro-ph/0012546) 
Zhang, Q., & Fall, S. M. 1999, ApJ, 527, L81 



